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O ■ Abstract 

'^'h ' We present the first data structures that maintain near optimal maximum cardinality and maximum 

f^ . weighted matchings on sparse graphs in sublinear time per update. Our main result is a data structure 

that maintains a (1 + e) approximation of maximum matching under edge insertions/deletions in worst 
case 0{^/m€~^) time per update. This improves the 3/2 approximation given in [Neiman, Solomon, 

r/^' STOC 2013] which runs in similar time. The result is based on two ideas. The first is to re-run a 

static algorithm after a chosen number of updates to ensure approximation guarantees. The second is 
to judiciously trim the graph to a smaller equivalent one whenever possible. 

We also study extensions of our approach to the weighted setting, and combine it with known 
frameworks to obtain arbitrary approximation ratios. For a constant e and for graphs with edge weights 
between 1 and iV, we design an algorithm that maintains an (1 + e)-approximate maximum weighted 
matching in 0{^/rn\og N) time per update. The only previous result for maintaining weighted matchings 

QQ , on dynamic graphs has an approximation ratio of 4.9108, and was shown in [Anand, Baswana, Gupta, 

(^ ■ Sen, FSTTCS 2012, arXiv 2012] 

m ■ 

p 

T::f 1 Introduction 

^^ ', The problem of computing maximum or near-maximum matchings in a graph has played a central role 
in the study of combinatorial optimization |LP86t IPS8 2]. A matching is a set of vertex-disjoint edges in a 
graph, and two variants of the problem are finding the maximum cardinality matching in an unweighted 

^ ', graph, and finding the matching of maximum weight in a weighted graph. The problem is appealing for 
several reasons: it has a simple description; matchings sometimes need to be improved by highly non-local 
steps; and certifying the optimality of a matching yields a surprising amount of structural information 
about a graph. On static graphs, the current best algorithms for maximum cardinality matching run in 
0{m^/n) time, on bipartite graph by Hopcroft and Karp |HK71| . and on general graph by Micali and 
Vazirani |MV80] . In the weighted case, algorithms with similar running times were given by Gabow and 
Tarjan [(ITQlj . and by Duan et al. [DPS 11] , 

A natural question from a data structure perspective is whether on a dynamically changing graph the 
solution to an optimization problem can be maintained faster than recomputing it from scratch after each 
update. For maximum cardinality matching, an 0{rn) time algorithm follows by executing one phase of 
the static algorithm described by Tarjan [Tar83j . For dense graphs, a faster running time of 0(n^'^^^) has 
been shown by Sankowski |San07] . and to date this is the only known result that gives sublinear time per 
update. For trees, Gupta and Sharma [GS09j gave an algorithm based on top trees that takes O(logn) time 
per update. 

On static graphs, a nearly-optimal matching can be computed much faster than finding the optimum 
matching. So it stands to reason that the same should apply in the dynamic case. Ivkovic and Llyod |IL94] 
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gave the first result in this direction: an algorithm that maintains a maximal matching with 0{{n + 
m)^ ) update time. Recently there has been a growing interest in designing efficient dynamic algorithms 
for approximate matching. Onak and Rubinfeld designed a randomized algorithm that maintains a c- 
approximation of maximum matching in 0(log n) update time |OR10j . where c is a large unspecified 
constant. Baswana, Gupta and Sen |BGSll| showed that maximal matching, which is a 2-approximation of 
maximum matching, can be maintained in a dynamic graph in amortized O(logn) update time with high 
probability. Subsequently, Anand et al. |ABGS12al IABGS12b] extended this work to the weighted case, 
and showed how to maintain a matching with weight that is expected to be at least 1/4.9108 « 0.2036 of 
the optimum. 

These results show that a large matching can be maintained very efficiently in dynamic graphs, but 
leave open the question of maintaining a matching closer to the optimum matching. Recently, Neiman 
and Solomon |NS12j showed that a matching of size at least 2/3 of the size of optimum matching can be 
maintained in 0{^/m) time per update in general graphs , as well as 0(logn/loglogn) time per update 
on bounded arboricity graphs. A similar result of maintaining 3/2-approximate matchings was obtained 
independently by Anand |Anal2| . This leads to the following question: Can we maintain a matching close 
to maximum matching (say (1 + e)-approximate matching) in a dynamic weighted or unweighted graph? 
We answer this question in affirmative by designing the first data structure that maintains arbitrary quality 
approximate max-cardinality and max-weighted matching in sublinear time on sparse graph. 

Our algorithm differs significantly from previous ones in that we do not maintain strict invariants. 
Baswana et al. [BGSllj maintained a maximal matching, which ensures no edge has both endpoints 
unmatched; and the 3/2-approximate algorithm designed by Neiman and Solomon [NS 12 j remove all length 
three augmenting paths in the graph at each update step. Our approach makes crucial use of the fact 
that the optimization objectives involving matching is stable. That is, a single update can only change the 
value of the optimum matching by 1. So if we find a matching close to maximum matching at some update 
step, it remains close to maximum even after several updates to the graph. In case the current matching 
ceases to be a good approximation of the maximum matching, we then re-run the static algorithm to get 
a matching that is close to optimum. This approach of re-running a expensive routine occasionally is a 
common technique in dynamic graph data structures |HK99l iHdLTOH IBKS12J . It is particularly powerful 
for approximating matchings since the stability property gives us freedom in choosing when to re-run the 
static algorithms. But re-running static algorithm occasionally works well when the maximum matching 
in the graph is large. To deal with graphs having small maximum matching, we introduce the concept of 
core subgraph which is the central concept of our paper. A core subgraph is a subgraph of a graph having 
the following two properties: Its size is considerably smaller than the entire graph. Secondly, the size of 
maximum matching in core subgraph is same as the size of maximum matching in the entire graph. We will 
crucially use these two properties in designing a dynamic algorithm for approximate matching. A detailed 
description of our algorithm, as well as other components of our data structure are presented in Section [3] 
and Appendix lAl The main result for approximating the maximum cardinality matching can be stated as 
follows: 

Theorem 1.1. For any constant e < 1/2, there exists an algorithm which maintains a (l + e)- approximate 
matching in an unweighted dynamic graph in 0{^/rn) worst case update time. 

It can be argued that the stability property of matchings that we rely on is rare among optimization 
problems. For most other problems like shortest paths and minimum spanning tree, there exist updates 
that require immediate changes in the approximate solution maintained. For matchings, such updates 
exists in the weighted version, where the objective is the sum of weight over edges in the matching. Direct 
extensions of our approach have linear dependencies on N in update time, where A^ is the maximum 
weight of an edge. This dependency can in fact be viewed as a quantitative measurement of the decrease 
in stability as we allow larger weights. 



As a result, we investigate rounding/bucketing based approaches which have logarithmic dependency 
on N in Section [H This was first studied for maintaining dynamic matchings by Anand et al. |ABGS12a] . 
and they used dynamic maximal matchings as a subroutine in their algorithm. Directly substituting our 
result for maximum cardinality matching leads to immediate improvements in the approximation ratio 
which is the second result in this paper. 

Theorem 1.2. For any constant e < 1/2, there exists an algorithm that maintains (3 + e)- approximate 
maximum weighted matching in a graph where edges have weights between [1,N] in 0{y/rnlogN) worst 
case update time. 

Our (3 + e)-approximation algorithm is derived from known schemes which bucket edges based on their 
weights. The rounding scheme we use in this algorithm is based on algorithm designed by Anand et al. 
[ABGS12a] . It is not clear whether any extension of this bucketing scheme will lead to a (l + e)-approximate 
matching. To do this, we devise a new rounding scheme which obtains arbitrarily good approximations of 
maximum matching, albeit at the cost of a much higher dependency on 1/e in the running time. 

Theorem 1.3. For any constant e < 1/2, there exists an algorithm that maintains (1 + e)- approximate 
maximum weighted matching in a graph where edges have weights between [1,N] in 0{y/m\ogN) worst 
case update time. 

As with the algorithm by Neiman and Solomon [NS12J , our algorithms are deterministic and the update 
time guaranteed by them is worst case. However, for simplicity in our presentations we will often start by 
describing the simpler amortized variants. 

2 Preliminaries 

We start by stating the notations that we will use, and reviewing some well-known results on matchings. 
An undirected graph is represented by G = {V, E), where V represents the set of vertices and E represents 
the set of edges in the graph. We will use n to denote the number of vertices \V\, and m to denote the 
number of edges \E\. 

A matching in a graph is a set of independent edges in the graph. Specifically, a subset of edges, 
M C £■ is a matching if no vertex of the graph is incident on more than one edge in M. A vertex is called 
unmatched if it is not incident on any edge in M, otherwise it is matched. Similarly, an edge is called 
matched if it is in M or free otherwise. A vertex cover is a set of vertices in a graph such that each edge 
has at least one of its endpoint in the vertex cover. 

The maximum cardinality matching(MCM) in a graph is the matching of maximum size. Similarly, 
given a set of weights w : E ^ [1,N], we can denote the weight of a matching M as w{M) = X^gg^ w{e). 
The maximum weight matching(MWM) in a graph is in turn the matching of maximum weight. We will 
use M to denote a optimum matching for either of these two objectives depending on context. 

For measuring the quality of approximate matching, we will use the notation of a-approximation, 
which indicates that the objective (either cardinality or weight) given by the current solution is at least 
1/a of the optimum. Specifically, a matching M is called a-MCM if \M\ > ^\MCM\, and a-MWM if 
w{M) > ^\MWM\. 

Finding or approximating MCMs and MWMs in the static setting have been intensely studied. Nearly 
linear time algorithms have been developed for finding (1 + e) approximations, and we will make crucial 
use of these algorithms in our data structure. For maximum cardinality matching, such an algorithm for 
bipartite graph was introduced by Hopcroft and Karp [HK71j . and extended to general graphs by Micali 
and Vazirani pVSOl IVazl2j . 



Lemma 2.1. There exists an algorithm ApproxMCM that when given a graph G with m edges along 
with a parameter e < 1, return an (1 + e)-MCM in 0{me~^) time. 

For approximate MWM, there has been some recent progress. Duan et al. |DPl(HlDPSlI] designed an 
algorithm that find a (1 + e) approximate maximum weighted matching in 0{me~^ log(e~^)) time. 

Lemma 2.2. IDPIC^ \DPSllf There exists an algorithm ApproxMWM that when given a graph G with 
m edges along with a parameter e < 1, return an (1 + e)-MWM in 0{me~^ log(e~^)) time. 

An logarithms in this paper are with base 2 unless mentioned otherwise. 

3 (1 + e)-MCMs Using Lazy Updates 

3.1 Overview 

To maintain approximate matching, we exploit the stability of the matching and use the static algorithm 
for matching AppROxMCM periodically. Our starting point is the observation that the size of maximum 
matching changes by at most 1 per update. This means that if we have a large matching that's close to 
the maximum, it will remain close to maximum matching over a large number of updates. So we use the 
following approach: Find a matching at certain update step and wait for certain number of updates till 
the matching is a good approximation of maximum matching. This approach works well if the maximum 
matching is itself large to begin with. But if the maximum matching itself is small, we still need to run 
the static algorithm many times. 

To overcome this, we show that instead of finding a maximum matching on the entire graph, we can 
use a small special subgraph such that the size of maximum matching in this subgraph is same as the size 
of maximum matching in the entire graph. We call this subgraph a core subgraph, and it is the central 
idea of our (1 + e) approximate algorithm. As this subgraph is considerably smaller, the time needed to 
find a maximum matching on it is considerably less. We will show that this core subgraph can be formed 
using the vertex cover of the entire graph. Specifically, we take the vertex-induced subgraph formed by 
the cover, along with some special chosen edges out of vertices belonging to the cover. 

But this leads to another question: How do we maintain a vertex cover in a dynamic graph? For 
this, we can use the algorithm of Neiman and Solomon |NS12] . One of the invariants in this algorithm 
is that there are no edges between unmatched vertices, which means the set of matched vertices form a 
2-approximate minimum vertex cover. Therefore reporting these vertices suffices for a vertex cover at any 
update step. However, note that our dependence on the above algorithm is not critical. Specifically, we 
design another simple algorithm which does not depend on the algorithm of Neiman and Solomon |NS12) 
for finding the core subgraph. A description of this, as well as modifications for handling edges with weights 
in a small range, and obtaining worst case bounds are in Section [3.31 with details deferred to Appendix lAl 

3.2 Algorithm 

We start with some notations that we will use in this section. We number the updates from 1 to t and use 
the following notations: 

• G{i): The graph after the i**^ update. 

• M{i): A matching computed on G{i) 

• M{i \j): Let delMii,j) denote the set of all edges in M{i) that are deleted from the graph between 
update steps i and j. We define M{i\j) to be M{i) \delM{i,j), i-e., M{i\j) consists of all the edges 
in the matching M(i) that are not deleted between update step i and j. 



Also, we will use Ai{i) to denote the optimal matching at step i. The approximation guarantees of 
M(i \ j) is as follows: 

Lemma 3.1. // e, e' < 1/2 and M{i) is an (1 + e)-MCM in G{i), then for j <i + e'|M(i)| M{i \ j) is an 
(1 + 2e + 2e')-MCM in G{j) 

Proof. Suppose there were /c^s insertions and k^ei deletions in the k = e'|M(i)| updates between updates 
i and j. The assumption about M{i) implies that |A^(i)| < (1 + e)|M(i)|. Since each insert can increase 
the size of the maximum matching by 1, we have |A^(j)| < |A^(i)| + kins- Also, each deletion can remove 
at most one edge from M{i), so \M{i\j)\ > \M\ — kdei- The approximation ratio is then at most: 

\M{j)\ ^{l + €)\M{i)\+hns 



\M{i\j)\- \M{i)\-kdel 

e\M{i)\+k 



1 + 



\M{i)\-kdel 



l/2\M(i)\ ^vacekdei<e\M{i)\<l/2\M{i) 

< 1 + 2e + 2e' 
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This fact has immediate algorithmic consequences for situations where the maximum matching is large. 
Suppose we computed an (1 + e/4)-MCM for G{i), M{i), then M{i \j) is (1 + e) approximate maximum 
matching as long as j < i+e\M{i)\/A. The 0{me~^) cost of the call to ApproxMCM (given by Lemma [2TT]) 
can then be charged to the next e|M(i)|/4 updates, giving 0{ ,^^-^, e~'^) time per update. When |M(i)| is 
large, this cost is fairly small. On the other hand, when |M(z)| is of constant size, this approach will make 
a call to ApproxMCM almost every update. 

For small size matching, we introduce the concept of core subgraph. As mentioned previously, core 
subgraph can be found by using a vertex cover G. 

Definition 3.2. Given a graph G and a vertex cover Vcover, o- core subgraph G' consists of: 

• All edges between vertices in Vcover 

• For each vertex v G Vcover , the \Vcover\ + 1 edges of maximum weight of v to vertices in V\ Vcover- In 
case of an unweighted graph, these edges can be chosen arbitrarily. 

An illustration of a core subgraph is shown in Figure [TJ It can be used algorithmically as follows. 

Lemma 3.3. Let G' be a core subgraph of G formed using a vertex cover Vcover ^ V . If M' is a {l+e)-MCM 
in G', then it's also a (1 + e)-MCM in G. 

Proof. We first show that the size of the maximum matching in G is the same as the size of the maximum 
matching in G' . Among all maximum matchings in G, let M be one that uses the maximum number 
of edges in E{G'). For the sake of contradiction, suppose M contains an edge {u,v) in E{G) \ E{G'). 
Since Vcover IS a vertex cover, one of u or f is in Vcover-, without loss of generality assume it's u. By the 
construction rule, for (u, f) to not be included in G", there exists |V"coi)er| + 1 neighbors oi u mV \ Vcover 
that are in G', let them be Ny\y^^^^^[u). As the maximum matching in G has size at most | Vcover I and 
there are no edges with both endpoints in y \ Vcover-, at most | Vcover I vertices in Ny\y^^^^^{u) can be 
matched. Therefore there exists an unmatched vertex x in N-y\y^^^^^{u) . Substituting (u, f) with {u,x) 
gives a maximum matching that uses one more edge in the G', giving a contradiction. 

Combining this with the fact that E{G') C E{G) implies that the size of the maximum matchings in 
G and G' are the same. Therefore any (1 + e)-MCM in G' is also a (1 + e)-MCM in G. D 




Figure 1: An example showing the core subgraph G' of G. All the vertices in the inner circle form a vertex 
cover Vcover of G. The core subgraph contains all the edges induced by the vertices in Vcover plus atmost 
|^oi;er| + 1 edges from each vertex v £ Vcover whose other endpoint is not in the vertex cover 



As mentioned previously, we can find Vcover in the graph by using the algorithm of Neiman and Solomon 
|NS12j . Their algorithm maintains 3/2 approximate matching in 0{^/rn) update time in the worst case 
which is less than the bound we are claiming. Whenever we need a vertex cover, we can report all the 
matched vertices in the 3/2 approximate matching. From now on we will assume an oracle access to the 
vertex cover at any update step. A more detailed treatment of maintaining a small cover can be found in 
Appendix lA.il 

Any vertex cover Vcover in graph G{i) formed out of a valid matching has the following property: 
l^overl ^ 2|A^(i)|. This is because the size of any valid matching is always less than the maximum 
matching size |A^(i)|. Therefore when |A^(i)| is small, we only need to run the static algorithm given by 
Lemma [2TT] on a core subgraph G'{i) of G{i). We can construct this graph in 0{\Vcover\^){= 0{\M{i)\^)) 
time by examining up to 0(|V^ot)er|) neighbors of each vertex in Vcover- Using Lemma l2.H we can find a 
(1 + e) approximate matching in this graph in 0(|A^(z)pe~^) time. Furthermore, Lemma |3. II allows us to 
charge this 0{\Ai{i)\'^ e~^) time to the next e|A^(z)|/4 updates. Therefore, cost charged per update can 



be bounded by 0{\M{i)\e'' 
following global states: 



'), which is small for small values of |A^(z)|. Our data structure maintains the 



1. A matching M. 

2. A counter t indicating the number of updates until we make the next call to ApproxMCM 

3. A vertex cover Vcover (Using the algorithm of Neiman and Solomon |NS12] ) 

Upon initialization M is obtained by running the static algorithm on G, or can be empty if G starts 
empty, t can be initialized to e/4|M|. Since we handle insertions and deletions in almost symmetrical 
ways, we present them as a single routine Update, shown in Figure [2] 

The bounds of this routine is as follows: 

Theorem 3.4. The matching M is an [l + e)-MCM over all updates. Furthermore, the amortized cost per 
update is 0{^/rne^'^). 

Proof. Let the current update be at time j, and the matching M that we maintained was computed in 
iteration i < j. So at update step i, the matching M is same as M(i) and at update step j, it is M{i \ j). 



1 if Update is a deletion and (u, v) £ M then 

2 Remove (u, v) from M; 

3 t ^ t — 1 ; 

4 if t < then 

5 Construct a core subgraph G' of the current graph; 

6 M ^ ApproxMCM(G', 1 + e/4) ; 

7 t ^ e/4|M| ; 



Procedure Update(u, v) 
Figure 2: Lazy update algorithm for maintaining (1 + e)-MCMs 

If t > 0, then since t was initiahzed to e/4|M(i)|, we have j — i < e/4|M(i)|. The guarantees for M{i \ j) 
follows from Lemma l3. II with e <— e/4 and e' ^ e/4. 

We now turn our attention to running time. Consider a call to ApproxMCM made at update i. 
Assume that e|M(i)| > 1. We have seen that there exists a core subgraph G'{i) such that the number of 
edges \E{G'{i))\ can be bounded by 0(min{m, |A^(i)p}). Since M(i) is a (l + e/4) approximate matching, 
(1 + e/4)|M(i)| > \M{i)\. So, the size of E{G'{i)) is 0(min{m, |M(i)|2}). Moreover, the cost of finding 
the matching( in ApPROxMCM) in the graph can be at most 0(min{?n, |M(i)p}e~^). This cost can be 
charged to the e|M(i)|/4 updates starting at update i, implying the following amortized cost per update: 

0(min{m, |M(i)|2}e-i) ^/ . f m ,,^,.^, 

If |M(i)| > -y/m, the first term inside min is at most ^/m, otherwise the second is at most ^/^m. Combining 
these two cases gives our desired bound. 

Now we take a look at some corner cases to complete the proof. We assumed that the cost of finding 
the matching at level i can be charged to next e|iVf(i)|/4 updates. This is true except for last call to 
ApproxMCM. The number of updates after this last call can be less than e|M(i)|/4. This cost can be 
amortized to all the updates. Since the number of updates is at least m, the total cost charged to each 
update step is 0(e~^). 

The other case is when e|M(i)| < 1. This implies that G'{i) has size at most 0(e~^) and finding a 
matching in such a graph takes time 0(e~^). For any constant e, this bound is 0{^/me~^) and can be 
charged to update step itself. 

So the amortized cost charged to any update step is at most 0{^/rne~'^). D 

3.3 Improvements, Worst-Case Bound, and Weights 

Several improvements can be made to the simpler version of our algorithm described above. We state the 
main statements here, and more details on these modifications can be found in Appendix lAl 

First, note that we depend on the algorithm of Neiman and Solomon |NS12| to maintain approximate 
vertex cover. Instead of using their algorithm, we design another simple dynamic algorithm which maintains 
approximate vertex cover. This algorithm is similar in spirit as our approximate matching algorithm, i.e, 
we use the property that vertex cover are also stable and a single update to the graph can change the 
vertex cover by 1. Using the techniques similar to the one presented in the previous section, we design 
an algorithm in Appendix I A. II which take 0{^/rn) update time in the worst case to maintain approximate 
vertex cover. 



Note that ApproxMCM may take 0{me~^) time in the worst case. So our algorithm in the pre- 
vious section had an amortized running time of 0(-\/me~^) per update. We show that we can maintain 
approximate matching in worst case 0{^/rne~'^) update time. Specifically, we show in Appendix IA.3I that 
computation cost of 0{me~^) time in ApproxMCM can be distributed across a number of updates. 

Furthermore, our ideas of maximum cardinality matching can also be adapted to maximum weighted 
matchings. This extension is natural because maximum cardinality matchings is a special case where all 
edges have weight 1 . A closer examination of the proofs of Lemma 13.11 shows that when all edge weights 
are in the range [1, N], the stability properties only degrade by a factor of A^. In Appendix lA.21 we present 
the following result: 

Theorem 3.5. For any constant e, there exists an algorithm that maintains (1 + e)- approximate maximum 
weighted matching in a graph where edges have weights between [1,N] in 0{^/'mNe~'^) update time. 

4 Approximate Weighted Matchings with Polylog Dependency on N 

We now show algorithms that approximate the maximum weighted matching in time that depends on 
log A instead of poly (A). This reduced dependency on A is a subject of study in static algorithms since 
A is often poly(n) or larger. 

Our overall scheme is based on the data structure for weighted matchings by Anand et al. |ABGS12al 
IABGS12b] . Their algorithm maintains log A levels and the edges are partitioned across various levels 
according to their weights. A matching Mi is maintained at each level /, and they gave a way to form a 
single matching M from these log A matchings. Algorithmically M can be viewed as the result of a greedy 
process which proceed in decreasing order of levels and adds edges whenever possible. Alternatively, it can 
be viewed as adding an edge {u, v) G Mi to M and removing all edges incident to u and v from all M//s 
where /' < /. At any update step, the matching maintained is equivalent to the one generated in Figure [3j 



M = 0; 


Let /max and /min be the maximum and minimum level number respectively; 


for I = /max to Irain do 




M = M\JMi; 




for (n, v) £ Mi do 




|_ Remove all the edges adjacent to u and v from Mi' such that /' < / 



Figure 3: Generating M 

Anand et al. |ABGS12al IABGS12b] showed that the combined matching M can be maintained on a 
dynamic graph if the matching at each level I can be maintained. We will use their result as a black-box 
via. the following Lemma. 

Lemma 4.1. (\ABGS12b\l ) If the matching on each level is maintained in 0{f{n)) update time, then the 
overall matching can be maintained in 0(/(n) log A) update time. 

In their work, /(n) = O(logn) due to the use of the dynamic maximal matching data structure by 
Baswana et al.[BGSTT], which leads to a total bound of 0(lognlog A). We will substitute our algorithms 
in place of this algorithm, and investigate different leveling schemes which lead to improved approximation 
ratios. This comes at a cost of a higher value of f{n) = O(i/mpoly(e^^)), which leads to a time of 
0(^/m log Apoly(e^-'^)) per update. 



In Section [^m we present a deterministic algorithm which maintains a (3+e)-MWM in 0{^/mlogNe~^) 
time, and in Section [4.21 we given an alternate approach which maintains a (1 + e)-approximate MWM 
in 0{^/mlogNe~'^'~'-^'^'' ') time per update. Note that in both the above algorithm, we will maintain 
approximate MCM or MWM matching at each level. For this we can use the amortized and worst-case 
versions of our data structures described in Sections [3] and Appendix |A] leading to corresponding types of 
final bounds for the above algorithm. 

In many of our proofs, we will incur (1 + 0(e)) multiplicative error in several places. As a result, the 
final approximation factors in our calculations will often be 1 + ce for some constant c. Such bounds can 
be converted to 1 + e approximations by initiating the calls with smaller values of e. As a result, we will 
omit these steps to simplify presentation. 

4.1 (3 + e)- Approximation Using Approx MCMs 

We first show that our data structure for maintaining (1 + e)-MCMs given in Theorem 11.11 can be used 
on each level. The transformation for turning a MWM problem into a set of 0(log A^) MCM instances is 
based on a rounding scheme by Eppstein et al. [ELMSlOl IELSW12] . For a fixed value of r, we assign an 
edge e with w{e) € [a'"*"^, a'"^''"^"'^) to level / where a is a constant which we will calculate later. Note that 
the level of some edges can be —1, but our proof below can extend to any negative level as well. We define 
the rounded weight of an edge e assigned to level / using: 

/ X def /_i_j. 

Wr[e) := a^ 

Our analysis of the quality of M is based on mapping each edge in M to a set of edges in M^'s. For 
e{u,v) G M from level /, we define TZ{e) as: 

7^(e) = {e} U {(x, y) \ {x, y) G Mi> where /' < /, and {x, y} D {u, v} / 0} 

In other words, TZ{e) contains edge e and all those edges adjacent to u and v from lower levels that were 
removed when {u, v) was added to M. Note that e is the only edge in TZ{e) from level i. And for all /' < /, 
there can be at most 2 edges from level /' in TZ{e). To simplify our notations, we will use tc(S') to denote 
the total weight of a set of edges S (that could be either M, Ai or Mi for some /') 

For an edge e G M, let $(e) denote the total rounded weights of edges in TZ{e), i.e., ^(e) = WriTZie)). 
We can show that $(e) is closely related to Wr{e). 



Lemma 4.2. For e G M , 



, , , a + 1 , , 

$ e < -Wr{e) 

a — I 



Proof. Let e G M be on level i. Since there are at most 1 edge on level i assigned to e (e itself) and 2 



edges per level assigned to e for each level j < i, we have: 



$(e)= Y. Me') 

e'e-R.(e) 

=Wr{e)+'Y ^ Wr{e') 

j<i e'eMjSze'enie) 

j<i 



<a*+" 


l + 2^a^-* 




< j<j 


=Wr{e) 


/ 1 \ 
1 1 'I 


r 1 ^ 1 


a + 1 
1 


u;r(e) 



D 
This allows us to relate the weight of M to the weight of the optimum matching, M . 
Lemma 4.3. 

{I + e)^^w{M) > Wr{M) 
a — 1 

Proof. Let M.{i) denote the edges of jH at level i. Since Mi is a (1 + e) approximate matching at level i, 
we have: 

\M{i)\ <{l + e)\Mi\ 
Wr{J^{i)) <(1 + e)'Wr{Mi) Since edges on same level have the same values of Wrie) 

Wr{M) <{l + e)YwriMi) 

i 

Consider an edge e = {u,v) G Mj. If e G M, then e E TZ{e). If e ^ M, then there exists an edge 
e' G M at level j > i such that one of the endpoints of e' is either u or v, which means e is in the set TZ{e'). 
Therefore each edge e can be mapped to one or more TZ{e'), and we have: 

^(M) >Ywr{Mi) 
i 

Which implies (1 + e)$(M) > Wr{M)- Summing Lemma 14.31 over all edges in M then gives: 

{I + e)^^Wr{M) > Wr{M) 

a — 1 
And the result follows from the fact that the rounded down edge weights satisfy Wr{e) < w{e). D 

Hence, it suffices to find ratio between Wr{M) and w{Ai). The analysis in Anand et al. |ABGS12b] 

bounded this ratio over a uniformly random choices of r. They showed that the expected rounded value 

a- 1 
of the optimum matching, Er[tfr-(A^)] satisfies Er[tUr(A^)] > — -. w{Ai), which when combined with 

a In a 
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Lemma [4. 31 leads to an expected approximation ratio of about 3 + e when a ~ 5.704. Here we show instead 
that a deterministic and worst-case bound can be obtained by using 0(l/e) versions of our data structure, 
each with a pre-selected value of r. 

We have k = lna/ln(l + e) copies of our algorithm which work exactly identically but with different 
value of r. For the j*^ copy, r{j) = ^^. Consider an edge e such that w{e) = a*"'" where < 6 < 1. Let 
j* is the value such that ^-^— < 6 < ^. Then we have: 

, , f 0'+'^ ifi<r 

Informally, an edge e is at level i in j*'^ copy, if j < j* otherwise it is at level i — 1. We want to relate the 
weight of maximum matching A^ in G to the new weight in these k copies. Specifically, we want to get 
a relation similar to the relation between Er[u;r(A1)] and w{M) mentioned above. We show that there 
exists a j with the following relation. 

Lemma 4.4. There exists a j such that: 

" ' a in Q 

Proof. Summing over all j of Wr(j){e) gives: 

1^ . ^^h ij-Lil 



EU^ride) _EU"'^^ +^J=,' 



+1 






w{e) «*+'' 

(j* k 

Since J* was chosen such that f > 5, -k6 + j*-l > -k{^)+j*-l = -land (l + e)-^'^+-'*-^ > {l + e)'^. 
Substituting this gives: 



> 



1 + ^rM E(i + ^)'"'" + E (1+^r'" 



The two summations is a rearranged version of a geometric sum. It can be rearranged by substituting 
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/ = J* — j + 1 and l=j*— j + k + 1 in the first and second summation respectively to obtain: 



■w{e) 



1=1 

l-(l + e)-'= 


(1 - 1/a) 


e 

a-1 



;=i*+i 



ae 
Summing this over all edges in M gives: 



eSAI j eeM 



By an averaging argument we get: 



k 
J 

aek 

Note that k = lna/ln(l + e). Here we make use of the following known fact about the behavior of the log 
function around 1: 

Fact 4.5. For e < 1, ifO<x<e, then ln(l + x) > (1 - e)x. 

Applying it with x = e gives: 

r /A^M (a-l)ln(l + e) 

m8iX.{Wr(i)(M)} =- -—^ -w(M) 

j " ' aelna 

^ {a-l){l-e)e ByFact|13] 

aelna 

a ma 

D 

Combining Lemmas 14.31 and 14.41 gives the following theorem. 

Theorem 4.6. For any e < 1/2, there exists a fully dynamic algorithm that maintains a (3 + e)-MWM 
for any graph on n in worst case 0{^/mlog Ne^"^) time per update. 
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Proof. Consider maintaining k copies of our data structure and picking the maximum weighted matching 
among these copies as the current best matching. 
Using Lemma 14.31 we get: 

Vi (1 + e)^w{M{j)) >Wri,)iM) 
a — 1 ' 



Using Lemma l4.4t there exists a j' = argmaXj{it;(M(j))} such that 



a - 1 

Wr(ii\{M) > (1 - e)— ; w{M) 

" ' a In a 



Combining the above two equations we get: 

(1 + e)'^w{M{j')) >(1 - e)^w{M) 
a — I a In a 

Where one can check that ^^ < (H-4e) when e < 1/2. By a suitable choice of e, this factor of l+4e can 
be turned into 1 + e'. This imphes that the approximation ratio obtained by our algorithm is - — '^^ _i\2° "° - 
This term achieves its minimum value of w 3 + 3e when a ~ 5.704. Again this approximation ratio can be 
turned into 3 + e' by a suitable choice of e. 

For the update time, note that since a is a constant, k = 0(l/log(l + e)) = 0(l/e) copies of the 
structure are needed. In each such copy, a matching can be maintained in 0(y^ log A^e~^) update time. 
So matching in all the copies can be maintained in 0{^/^m\ogNe~^) time per update. D 

4.2 (1 + e)-MWMs Using Approximate MWMs 

Overview 

In this section, we present an algorithm that maintains a (1 + e)-MWM using a more gradual bucketing 
scheme. We start by observing the definition of lZ{e) for an edge e{u, v) in M from the previous section. 
Informally, TZ{e) contains edge e and all those edges adjacent to u and v from lower levels that were 
removed when (n, v) was added to M. A closer look at our algorithm reveals that the approximation ratio 
depends on the ratio of weight of e and the combined weight of edges in 7l{e). This ratio can be reduced 
if the edges at lower level have significantly less weight than the weight of edge e. To achieve this, we will 
artificially create levels such that the ratio of weight between two consecutive level is significant. For this, 
we will drop some edges from the graph to create a gap between two consecutive levels. In order to account 
for the weight of these dropped edges, we in turn need to keep several copies of our data structure with 
different edges left out in the other copies. 

We then proceed with the same algorithm as mentioned in Section U] with one main difference. Instead 
of maintaining (1 + e)-MCM at each level, we maintain (1 + e)-MWM at each level using the Theorem 13.51 
Note that this theorem has a dependence of A^ in its running time. We will show that each level can be 
formed in such a way that N can be bounded by 0{e~^^'' '). So the running time for maintaining (1 + e) 
approximate MWM at each level will have exponential dependence on (1/e). 

Thereafter, we combine the matching across the various level using the same procedure as mentioned in 
Section |H We will show that there exists a copy of our data structure such that the weight of the matching 
maintained by our algorithm in that copy is a good approximation of maximum weighted matching in the 
entire graph. 
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Figure 4: Bucketing and level scheme where C = 3 



Algorithm 



Once again we partition the edges by weights geometrically: an edge e is in bucket b if w^e) is in the 
range [e~ , e~^ ''~^^). However, our levels no longer corresponds to individual buckets, but instead to a set 
of C — 1 continuous buckets for value of C to be specified. We will also remove some of these buckets, and 
the choices of buckets to remove leads us to run several copies of our data structure simultaneously. 

We will run C = [e~^] copies of our algorithm, where in the c copy, we remove all buckets i such 
that i mod C = c. This leads to a set of graphs G^ . . .G'-'''^. Removing the buckets creates natural 
partitions of the remaining edges, which gives our levels. For a copy c, we will place buckets with b = 
[IC + c + 1 . . . (/ + 1)C + c — 1] into level I. Note that the ratio of maximum to minimum edge weight in 
each level is bounded by e^^'-^^^'{= 0{e^^^'' )). Therefore, the algorithm given in Theorem 13.51 allows 
us to maintain an (1 + e)-MWM in 0(e~*^(^ ^ ^/me~'^ log{e~^)) = 0{€~'^~'^^^~ ^ ^/rnlog{e~^)) time at each 
level. These matchings can in turn be combined together in the same way as in Section [H An illustration 
of levelling scheme used by our algorithm is shown in Figure HI 

We start by analyzing the guarantees of our algorithm on the c**^ copy. Specifically, the approximation 
ratio of the combined matching M^ w.r.t. the maximum matching Ai^ in this copy . Let M^, and Aif be 
the edges of M^ and Ai^ at level / respectively. Once again, for an edge e = {u,v) in Mf we define Tl'^ie) 
as: 

7^'=(e) = {e}U{(x,y)|(x,y) G Mp where /' < /, and {x,y}r\{u,v} / 0} 

For an edge e S M'^, let <l>'^(e) denote the total rounded weights of edges in TZ^{e), i.e., <l>^(e) = w{TZ^{e)). 
We can show that ^'^{e) is related to w{e) by the following inequality: 

Lemma 4.7. For any edge e in the combined matching M^, we have: 

^'(e) < (l + 3e)w(e) 
Proof. Assume that e is on level I. Since there are at most 1 edge on level I assigned to e (e itself) and 2 
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edges per level I' < I assigned to e, we have: 

e'e7^=(e) 

=wie)+Y, Yl "'(^') 

l'<l e'eM^fHWie) 

Since an edge on level I' is in bucket [I'C + c + 1 . . . {I' + 1)C + c — 1] and an edge in bucket b has weight 
at most (l/e)''"^^, the weight of an edge at level I' is < e""^"" +i)C_ 

i'<i 

Which can in turn be bounded relative to w{e). Since e is in level / and an edge in bucket b has weight at 
least (l/e)^ w{e) > g-'t^-^-i 

$'(e) <u;(e) + ^2e('-''~i)^+iu;(e) 
i'<i 

2e 

<w[e) + (jW{e) 

<u;(e)(l + 3e) assuming that e < 1/2 

D 
We now show the relation between the combined matching M'^ and Ai^. 
Lemma 4.8. 

(1 + 7e)w{M'') > wiM") 
Proof. Since w{M^) is a (1 + e) approximate maximum weighted matching on level /, we have: 

w{Mr)<{l + e)w{M[) 

w{M^) <(l + e)^u'(Mf) 
I 

Consider an edge e = {u,v) £ Mj^. If e G A4^, then e G 7lP{e). If e ^ IVP, then there exists an edge 
e' at level /' > / such that one of the endpoints of e' is either u or v, which means e is in the set TZ^{e'). 
Therefore each edge e can be mapped to one more TZ'^{e'), and we have: 

I 

Which implies (1 + e)^'^{M'^) > w{M^). Summing Lemma 14.71 over all edges in M'^ then gives: 

(1 + e)(l + 3e)w{M'') > wiM") 

And the bound fohows from (1 + e)(l + 3e) = 1 + 4e + 3e^ < 1 + 7e when e < 1. D 

We now find a relation between ^^{M^) and w{Ai). We show there is at least one copy whose maximum 
matching has weight at least (1 — 1/C) of w{Ai). 
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Lemma 4.9. At any update step, if the maximum weight m,atching in the current graph is w{A4), there 
exist a copy c such that w{M.'^) > (1 — l/C)w{M.). 

Proof. Let M'^ denote the set of edges in Ai that are not present in the c*'^ copy. Since each bucket is 
removed in only one copy, we have: 

iJcM" =M 
"^w{M^) =w{M) 

c 

Since Ai \ Ai'^ is a matching in C^, we have w{Ad^) > w{A4) — w{Ai^). Note that the inequahty is due to 
A^*^ being the maximum weighted matching in G^ instead of the restriction of A4 on it. Summing over ah 
c copies gives: 



c c 

=C-w{M)- ( J^t«(7W")] 



={C-l)-w{M) 

Dividing both sides by C gives that the average of w(A4^) is at least (1 — l/C)w{A4). Therefore there 
exist some c where w^Af^) > (1 — l/C)w{Ai). 

a 

Combining Lemmas 14.91 and 14.81 we deduce that there exists a copy c such that 
(1 + 7e)w(M'=) > (1 - e)w{M). So, 

w{M) ^ (1 + 7e) 
w(M^) ~ (1 - e) 

This ratio is less that (1 + 16e) for e < 1/2. By a suitable chose of e', the factor of (1 + 16e) can be turned 
into (1 + e'). 

This means that if we set C = \e~^~\ and maintain (1 + e)-MWMs on each copy of our data structure, 
then one of the maximum weight matching among these C copies will always be a good approximation of 
the maximum weighted matching for the entire graph. 

Note that in each copy of our data structure there are 0(log^-i N)/C = 0(^=q^%-rrr) levels and in 

each level an approximate MWM is maintained in 0{e~'^~^^^> ^/^m\og{e~'^)) time. This implies that the 
overall update time taken by our algorithm across the C copies is 0{C ■ cio^fe-^) \P^^~'^~^ lc)g(e~^)) = 
0{y/m,e~'^~^^'^> logA^). So we can state the following theorem: 

Theorem 4.10. For any e < 1/2, there exists a fully dynamic algorithm that maintains a (1 + e)-MWM 
in worst case 0{^/me~'^~^^'''' log N) time per update. 

5 Conclusion 

We showed a simpler method for maintaining approximate matchings that maintains (l+e)-approximations 
in about ^/rn time per update. Natural directions for future works are whether the update time can be 
improved, and whether the exponential dependency on e~^ in the weighted case can be reduced. Also, it 
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would be interesting to explore whether this type of approach can also maintain approximations of other 
optimization objectives. 

Theoretically our arbitrary quality approximation algorithm from Section 14.21 outperforms the (3 + e) 
approximation given in Section 14. 1[ However, its exponential dependency on e~^ makes it fall short of a 
practical algorithm for maintaining (1 + e)-MWMs. We believe that a more intricate rounding scheme such 
as the one given in Section 14. 2|, or possibly a data structure that incorporates details of the Duan et al. 
algorithm [DPIOI IDPSllj are promising approaches in this direction. 
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A Improvements, Worst-Case Bound, and Weights 

We now give the full details for the improvements described in Section 13.31 This includes removing the 
dependency on other data structures for finding Vcover in Section lA.il incorporating weights in Section P\,.2l 
and obtaining worst case guarantees in Section [A. 31 Each of these improvement can be made independent 
of the other ones. However, we will only give the pseudocode for a version incorporating all three extensions 
in Section IA.3I 

A.l Simpler Maintenance of Small Vertex Cover 

We now present a simple algorithm to maintain a constant approximation of vertex cover to remove our 
dependence on another algorithm in Section [3l Like our algorithm for maintaining approximate matching, 
a natural approach for maintaining Vcover is to use a static algorithm, ApproxCover for it and then 
do nothing for certain number of update step till the current vertex cover is a good approximation of 
minimum vertex cover. The AppROxCovER routine finds a vertex cover by finding a maximal matching 
and reporting all the endpoints of matched edges as the vertex cover of G. This gives a cover whose size is 
at most twice of optimum, and we can also bound its guarantees after a number of updates. We will use 
Vcoverii) to denote a vertex cover generated by calling ApproxCover(G(z)). Let M{i) be the maximal 
matching which is used to find the vertex cover Vcover (i)- In order to use this cover at some subsequent 
update J, we insert all vertices involved in edge insertions from update step i to j and denote this new 

cover as Vcover {i -^ j)- 

Lemma A.l. If Vcover{i) is a 2-approximate minimum vertex cover on G{i) and j <i + l/4M(i), then 
Vcover{i -^ j) is an 5 -approximation to the minimum vertex cover in G{j). 

Proof. We use that fact that vertex cover after t updates can decrease by at most t. The guarantee of 
ApproxCover implies that the size of the minimum vertex cover in G{i) is at least \Vcoverii)\/'^- Since 
M{i) < \Vcoverii)\, the size the minimum cover in G{j) is at least 

"^l ''cover \'^)\ "tI ''coDerv^JI — tI "^coiier (^^j| 

Also, since M{i) < \Vcover{i)\, 

\Vcover{i -^j)\< \Vcover{i)\ + l/4M(i) < 5 / 4\Vcover {i)\ 

Therefore \Vcoverii — ^ j)\ is a 5-approximation to the minimum vertex cover in G{j). D 

To maintain a small vertex cover, we find a new vertex cover using Vcover{i -^ j)- Similar to Lemma [3.3t 
we can build a core subgraph G' by taking a small set of neighbors of each vertex in Vcover{i -^ J)-, and 
prove the analogous result for maintaining vertex covers. 

Lemma A. 2. Let Vcover{i — ^ j) be a vertex cover in G{j) generated as above from Vcover{i), which was in 
turn returned by a call to AppROxCovER(G(i)). Consider the core subgraph G' having following edges. 

• All edges between vertices in Vcover{i -^ j) 

• Up to \Vcover{i " ^ i)l + l arbitrary edges to vertices in V\Vcover{i — ^ j) for each vertex in Vcover{i — ^ j)- 

ApproxCover(G') returns a new vertex cover for G{j), Vcoverij) whose size is at most twice the 
minimum. 
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Note that although Vcoverii -^ j) is a vertex cover in G{j), its may not be as good of an approximation 
due to adding all vertices involved in insertions between updates i and j. Generating a new cover can be 
viewed as a way to reduce the error. 

Proof. Let M' be a maximal matching found on G' by ApproxCover. It suffices to show that M' is 
a maximal matching on G{j) in a way similar to our proof of Lemma 13.31 For contradiction, suppose 
there exists such a free edge (u, v) in G{j) with respect to matching M' such that both u and v a are 
unmatched. As Vcoverii — ^ j) is a vertex cover for G{j), we may assume that u G Vcoverii -^ j) without 
loss of generality. The maximality of M' in G' implies that {u,v) ^ E{G'), and therefore v ^ V{G'), 
which means u has at least \Vcoverii — ^ j)| + 1 neighbors in F \ Vcoverii — ^ j)- As there are no edges 
between vertices in y \ Vcoverii -^ j), each edge in M' can only match one vertex from this set. Also, 
\M'\ < {Vcoverii -^ j)\ since the size of any matching is at most the size of a vertex cover. This means u 
has at least one unmatched neighbor in G' , contradicting the assumption that M' is a maximal matching 
in G'. a 

Using Lemma lA.ll and IA.21 we see that the vertex cover in graph G' can be found in Oi\Vcoverii -^ 
j)P)(= 0(|M(j)p)) time. This time can be amortized over the next M(j)/4 steps. Similar to the analysis 
in Section [3l we can show that the amortized update time at each step is Oiy/m) per update. 

A. 2 Incorporating Weights 

Till now we have described dynamic algorithms for maintaining approximate matching in unweighted 
graphs. Now we show that the same technique can be used to maintain maximum weighted matching. In 
order to handle weighted matchings, we will substitute the weighted matching algorithm given in Lemma[2?2] 
as our static algorithm. We also need to prove the analogs of Lemmas 13.11 and 13.31 for the weighted setting. 

Lemma A. 3. // all edge weights are in the range [1,N] and M is an (1 + e)-MWM in G'(i), then for 
j <i + j^\M\, M(i \ j) is an (1 + 2e + 2e')-MWM in Gj 

Proof. Similar to the proof of Lemma |3. 11 except each insertion can increase wiAi) by at most N, and 
each deletion can decrease wiM.) by at most N. D 

Lemma A. 4. Let Vcover ^V he a vertex cover, and G' be the core subgraph consisting of: 

• All edges between vertices in Vcover 

• The top I Vcover \ + 1 maximum weighted edges to vertices inV\ Vcover for each vertex in Vcover ■ 

If M' is ail + e)-MWM in G' , then it's also a (1 + e)-MWM in G. 

Proof. Similar to Lemma 13.31 but we need to compare the weights of edges (u, x) and (n, v) where x £ 
Nv\Vcover (^) ^'^'^ '^ ^ Vcover ■ By our choice of the | Vcover I edges incident to u being the ones of maximum 
weight, we have w(n, x) > wiu,v). D 

These two lemmas allows us to adapt the algorithm from the Section [3] for weighted matchings, with 
an additional factor of 0(-/V) in the running time. 

To construct G' , for each u S Vcover we need to examine edges incident to it in decreasing order of 
weights until we have either exhausted the list, or found |Vcoi,er| + 1 ones in 1/ \ Vcover- In either case, it 
suffices to extract the 2\Vcover\ + 1 edges of maximum weight incident to u. 
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A. 3 Worst Case Runtime 

In this section, we will present an algorithm which obtains a better worst case bound for maintaining 
approximate matchings. The algorithm is similar in spirit to the one described in Section [3l That is, we 
first find a approximate weighted matching in graph G{i) and then do nothing for the certain update step j 
till the matching we had obtained gives a good approximation of optimum. In order to find an approximate 
weighted matching in G(j), we perform the following tasks: 

• Finding a new vertex cover Vcoverij) using Vcover{i — ^ j) ( Lemma [A. 2p 

• Constructing core subgraph G' from VcoveriJ) (Lemma IA.4P 

• Running AppROxMWM(G', e) 

As G' may be as large as G, the number of edges in G" in worst case can be as high as 17 (m), and 
the call to ApproxMWM can only be bounded by 0{me~^ log(e~^)). However, Lemma I A. 3 1 implies that 
once an (1 + e/4)-MWM is found, it will remain an (1 + e)-MWM after e|M|/47V updates. Note that the 
algorithm does nothing in these e\M\/4:N update steps. In order to obtain a worst case bound, we perform 
useful computations in these update steps. Specifically, instead of finishing the three task mentioned above 
at one update step, we gradually work on these three task for these e\M\/AN steps. 

We now describe our algorithm in Figure [3 Our algorithm works in rounds. A round consists of three 
tasks that we mentioned above: 



1 if Update is a deletion then 

If (n, v) G M, remove it from M; 

Notify call to ApproxMWM in current round to remove (u, v) 

Notify current call to ApproxCover to add u. 

5 Perform 0{N ^/me~'^ log{e~^)) steps in the current round; 

6 if Current round finishes then 
Replace M by the results of the call to ApproxMWM; 
Initiate the next round consisting of: 



1. Construct core subgraph G' using the result from the previous call to ApproxCover 

2. ApproxMWM(G', 1 + e/8) 

3. ApproxCover(G') 



Procedure Updatelmp(u, v) 
Figure 5: Maintaining (1 + e)-MWMs with Worst Case Guarantees 



Constructing core subgraph G' from a vertex cover 

We obtain this vertex cover from the result of ApproxCover(G') which was called in the previous 

round. 

Running AppROxMWM(G', e) 

The result of this procedure is used as the approximate maximum weighted matching returned by 

the algorithm for the next round. 
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• Constructing a new vertex cover on G' using ApproxCover(G') 

The vertex cover returned by this procedure wih be used for the next round. 

To obtain a worst case bound, we will run these procedure gradually, i.e, we will do only 0{N^/rne~'^ log(e~^)) 
amount of work at each update step. However, this approach has a drawback. For example, suppose that 
while computing AppROxMWM(G', e), we added an edge e to the matching at some update step. But sub- 
sequently, if this edge is deleted from the graph, we don't have any way of reflecting this change in the future 
iterations of ApproxMWM(G', e). So we just remove these edges from the result of APPR0XMWM(G', e) 
(See line[3|). Similar is the case for the procedure ApproxCover(G'). However we will show that we don't 
lose much in terms of approximation due to this workaround. 

We will show that the number of update steps in a round are large enough so that all these three tasks 
can be performed if the work done at each update step is 0{N^/me~'^ log{e~^)). On the other hand, we 
also show that the number of update steps is small enough so that the approximate weighted matching 
returned in the previous round has a good approximate ratio till the current round ends. 

We find the the number of updates it takes for a for a round to finish using a proof analogous to that 
of Theorem 13.41 

Lemma A. 5. // a round initiated at update i with a vertex cover Vcover satisfying \Vcover\ ^ lQw{M{i)) , 
and is run for 0{N^/me^'^\og{e~^)) steps per subsequent update, it will finish within "^"^ gjv updates. 

Proof. The number of edges in G' can be bounded by min{r?T,, 0(|T4o»;erP)} = 0(min{77i, w{M.{i)) }). This 
means that the total cost of the current round is 0(min{r?T,, w(A^(i)) }e^^ log(e~^)). For this to finish in 
'^'"gjY updates, the amount of work that needs to be performed at each update is: 

O {rn:.n{mMM{i)f}e-'\og{e-')) / (^^^^^^ 



=0 I^Nm\n{-^j^^MM{i))}e~'\og{e- 

Considering cases of w{M.{i)) > \/m and tt;(A4(i)) < ^Jrn leads to an 0(N \/ral<?\og(er^y) worst case 
time per update step. D 

The worst case bounds can now be obtained by analyzing the duration for which M(i \ j) is used as 
M in the algorithm, and proving inductively the conditions required for the vertex cover in Lemma lA. 51 



We finish this section by giving the proof of Theorem 13.51 

Proof. We first show inductively that at the start of a round initiated at update j we have a vertex cover 
Vcover satisfying |yco»;er| < Ww{Ai{j)). The base case follows from both being 0, and for the inductive 
case, consider a round initiated at time j. Suppose the current value of Vcover was produced by a round 
initiated at time i < j. So the vertex cover Vcover at update step i is Vcover (i) and at update step j, it 
is Vcoverii " ^ j)- Lemma [A. 2 1 implies that Vcoverii) is a 2-approximate vertex cover on G{i). Using the 
conditions in Lemma f A. 5 1 we have j — i < ^""gTv < l/4|A^(i)|. Therefore by Lemmas lA.H |14oDer(^ — ^ j)\ 
is a 5-approximation of minimum vertex cover in G{j). Since the size of minimum vertex cover in G{j) is 
< 2|A^(j)|, this implies that {Vcoverii -^ j)\ < 10|A4(i)|. Thus the requirements of Lemma [A. 5 1 is satisfied 
for all rounds. 
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For the guarantees on M, it suffices to prove that the matching produced by a round initiated at update 
i remains an (1 + e)-MWM until the completion of the next round. Let the update by which this round 
finishes he i' . Applying Lemma lA.51 to all rounds gives: 

' - + 8iV 

and that the round initiated at update i' finishes in "^"^ gjv updates. Therefore it suffices to show that 
M{i \ j) is a (1 + e)-MWM for Gj for all j <i + ^{w{M{i)) + w{M{i'))). 

Since the call to ApproxMWM is ran with error e/4, we have w^M^i)) < (1 + e/4:)w{M{i)), and: 

w{M{i')) <w{M{i)) + N{i' - i) 
<{1 + €/8)wiMii)) 
<(l + e/8)(l + e/4)u;(M(i)) 

Summing these then gives: 

^{w{M{i)) + w{M{i'))) <^((1 + e/4)w{M{i)) + (1 + e/8)(l + e/A)w{M{i))) 

=^(2 + e/8)(l + 6/4KM(i)) 
<^wiM{{)) 

The bounds then follows from Lemma lA.31 with e' = 3e/8. D 

One other minor issue is that the next round need to access a previous state of the graph. As there 
is at most one round active at once, we can simply keep two versions, the current one as well as the one 
from when the round was initiated. A conceptually simpler way to address this issue can also be found in 
persistent data structures |DSST89] . which allows access to previous versions of the graph. 
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